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We report on the magnetic field (0T< B < 9T) dependence of tfie longitudinal thermal conduc- 
tivity K(r, B) of highly oriented pyrolytic graphite in the temperature range 5 K < T < 20 K for 
fields parallel to the c— axis. We show that k(T, B) shows large oscillations in the high-field region 
(B > 2 T) where clear signs of the Quantum-Hall effect are observed in the Hall resistance. With 
the measured longitudinal electrical resistivity we show that the Wiedemann-Franz law is violated 
in the high-field regime. 



INTRODUCTION 

Over the last 60 years the literature contains a substan- 
tial number of measurements and theoretical work on the 
electrical and thermal transport properties of the semi- 
metal graphite 1, 2]. In contrast to the common believe, 
however, several transport properties of graphite are not 
understood and some theoretical assumptions done in the 
past seem now less plausible. These doubts have their 
origin in the relatively new and controversial physics of 
a two-dimensional (2D) electron gas. Graphite has been 
tacitly assumed to be a 2D material, but the relatively 
low quality of the samples prevented clear measurements 
of its true 2D transport properties. The number of open 
questions regarding the transport properties of graphite 
is significant. Already early work noted that the mag- 
netic field dependence of the electrical longitudinal resis- 
tivity pxx and Hall resistance Rh^ even including ad- hoc 
dispersion relation and trigonal warping of the constant 
energy surfaces, was (and still remains) basically unex- 
plained Q- 

Recent measurements of the longitudinal resistivity of 
highly oriented pyrolytic graphite (HOPG) samples show 
a clear magnetic field driven metal-insulator transition 
(MIT) with a giant magnetoresistance at low fields and 
at low temperatures 0. For example, the longitudinal 
resistivity at 4.2 K can increase by more than one order 
of magnitude with a field B ~ 0.2 T. This MIT shows 
a scaling as found for 2D electron systems with similar 
scaling exponent but with a critical field ~ 0.1 T ap- 
plied perpendicular to the graphene layers Possi- 
ble origins for the MIT in graphite are being discussed 
nowadays in the literature in terms of superconducting 
fluctuations 0, excitonic insulator state triggered by a 
magnetic catalysis phenomenon and/or a Bose metal 
transition 6., _S] . High-resolution angle dependence of the 
magnetoresistance along the c— axis in HOPG samples in- 
dicates that the transport between layers gets incoherent 



the better the sample quality - characterized by the full 
width at half maximum (FWHM) of the rocking curve - 
and suggests that the coupling between planes is much 
less than the commonly assumed 0.3 eV [2|. 

In the quantum limit (QL) when only the lowest Lan- 
dau levels of graphite are occupied (B > 2 T), the lon- 
gitudinal electrical resistivity shows a reentrance to a 
metallic-like state below a field-dependent temperature 
T„r[B) mini. The function T„(B) depends on the di- 
mensionality of the graphite sample and it oscillates as 
a function of field for quasi 2D samples. This behavior 
might be an evidence for field-induced superconductiv- 
ity at the QL in graphite, discussed in Refs. [l2|,li3]- In 
the QL the Hall resistance Rh shows clear signs of the 
quantum Hall effect for samples with small FWHM [l]| . 

For ideal graphite the conduction electrons are ex- 
pected to follow a Dirac dispersion relation. These quasi- 
particles (QP) should have some similarities with the 
nodal QP of the d-density wave (DDW) state, which 
applies also to high-temperature superconductors. For 
a not superconducting DDW QP the violation of the 
Wiedemann-Franz (WF) law ,14j has been predicted 
as a characteristic of the relativistic spectrum of the QP. 
Briefiy, in the DDW state the QP can carry electrical cur- 
rent much more effectively than thermal current in the 
limit of very small QP density, a situation that applies 
to graphite (Fermi Energy Ep ~ 200 K). The theoreti- 
cal work of Ref. however, goes a step further taking 
into account the presence of a magnetic field within the 
DDW state in the calculations. This work finds that the 
WF law is restored at low enough temperatures (com- 
pared with Ep) independently of the applied field. We 
note that in none of these references quantum Hall states 
have been taken into account in the calculations. Taking 
into account recently published experimental evidence 
jlll |. these theories cannot be applied straightforwardly 
to graphite. The aim of our work is to check whether 
the WF law applies in graphite in particular in the QL 
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regime. 

The magnetic field dependence of the thermal conduc- 
tivity of strongly anisotropic HOPG samples was already 
measured in the past but one does not find curves in the 
published literature with the necessary resolution in the 
QL regime. There are at least two unpublished studies 
which deserve our attention. In his Ph.D. thesis, Ayache 
[l7t | measured k{T,B) and observed well defined oscilla- 
tions of k{B) at fields above 2 T and at constant tem- 
peratures below 10 K. This behavior was apparently also 
observed by WooUam jl^] but the curves were not in- 
cluded in the corresponding paper. Interestingly, both 
authors emphasized that the oscillations in k{B > 2 T) 
as a function of field were apparently in phase with the 
electrical resistivity oscillations, although the electronic 
contribution to k at high fields freezes out according to 
the Wiedemann-Franz law. In spite of this striking re- 
sult neither the curves nor any analysis of the data was 
published to our knowledge. 

High-resolution measurements of the field dependence 
of the thermal conductivity of quasi 2D HOPG samples 
have nowadays special relevance. In particular because 
evidence for a quantum-Hall behavior in 2D samples in 
the QL has been recently reported ^3] ■ The behavior of 
k(T, B) in graphite is not only important to understand 
the nature of the QP in this material at the QL but also 
gives us the chance to study the thermal transport in the 
quantum-Hall effect (QHE) regime. We note that the 
thermal transport in the QHE regime still remains an 
unclear problem since the discovery of this effect. In this 
work we have measured the longitudinal thermal con- 
ductivity K of a well characterized HOPG sample as a 
function of magnetic field applied parallel to the c— axis 
of the graphite structure. These measurements were ac- 
companied by measurements of the longitudinal and Hall 
resistances. The results clearly show that the WF-law is 
violated in the QL, whereas deviations are observed at 
lower fields. 



EXPERIMENTAL DETAILS 

The HOPG sample from the Union Carbide company 
measured in this work was selected because of its clear 
quasi-2D properties observed in the electrical and Hall 
resistivity. This behavior is correlated with the small 
FWHM= 0.26° of the rocking curve, which is a measure 
of the misorientation relative to the c— axis of the crys- 
tallites in the sample. The measured FWHM is one of 
the smallest we have obtained for HOPG samples. We 
note that most of the graphite samples studied in the 
literature of the 70's and 80's had much larger FWHM. 
Our experimental evidence indicates that in samples with 
FWHM values larger than ~ 0.5° the transport proper- 
ties show clear sign of 3D behavior and coherent trans- 
port in the c— direction]^ and therefore they do not 



reflect true 2D behavior of ideal graphite. The dimension 
of our sample was: length = 1.6 mm, width = 1.2 mm 
and thickness ~ 60 /im. The room temperature out-of- 
plane/basal-plane resistivity ratio at -B = of the sample 
was Pel Pb ~ 5 X Iff*. Furthermore, this sample does not 
show any maximum in the c— axis resistivity as a func- 
tion of angle for fields parallel to the graphene planes . 
This indicates incoherent electrical transport expected 
for weak-coupled conducting planes and for samples with 
a low density of defects. 

The Hall resistance was measured using the van der 
Pauw configuration with a cyclic transposition of cur- 
rent and voltage leads [ijl |20| at fixed applied-field po- 
larity as well as magnetic field reversal; no difference in 
Rh {H, T) obtained with these two methods was found. 
For the measurements, silver past electrodes were placed 
on the sample surface, while the resistivity values were 
obtained in a geometry with an uniform current distri- 
bution through the sample cross section. All resistance 
measurements were performed in the Ohmic regime. The 
absolute value of the longitudinal (basal-plane) resistiv- 
ity at zero field was pb(6 K, 0) ~ 2.4/ir2cm and 2.6/ir2cm 
at 10 K. The error in the absolute value is estimated to 
be ~ 30% due to geometrical errors. At 10 K pb reaches 
its low-temperature rest value within 10%. The resistiv- 
ity increases by two orders of magnitude for an applied 
field of 1 T at T < 10 K. 

For the longitudinal thermal conductivity measure- 
ment the temperature gradient (of the order of 200 to 
300 mK) was measured using a previously field- and 
temperature-calibrated type E thermocouples with a dc- 
picovoltmeter |2 j |. The thermocouple ends were posi- 
tioned one at the top and the other at the bottom of 
the main surface of the sample. A detailed calibration 
below 8 K was performed because in this temperature 
region the thermopower of our thermocouple is specially 
sensitive on the magnetic field with a non-monotonous 
dependence j2 ij . The experimental arrangement was re- 
cently used to study the longitudinal and Hall ther- 
mal conductivities of high-temperature superconducting 
crystals [23, [13 ■ We note that in general the measured 
thermal conductivity is k = — TS^cJi, where Ki is the 
"real" thermal conductivity of the sample, S the ther- 
mopower and Ui the electrical conductivity. In the case 
of our graphite sample the correction term to n is four 
orders of magnitude smaller than Ki at 10 K. 

Our system enables us to measure k{B) with a relative 
resolution better than 0.1% above 5 K. The thermal sta- 
bility was better than 10 mK in the whole temperature 
range 5 K < T < 20 K and magnetic field T< B < 9 T. 
The absolute error in the thermal conductivity was esti- 
mated to be < 30%. The obtained absolute value of k and 
its temperature dependence are similar to those from pre- 
vious studies |i UtI 111 El list . For example, at 10 K we 
obtain k(0) ~ 130 W/mK and k(0) ~ 33 W/mK at 5 K 
and zero fields. In all magnetic-field runs k{B) showed 
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reversible behavior. Irreversible behavior has been ob- 
served but we could prove that it was related to small 
temperature drifts since in the temperature range of the 
measurements k depends strongly on temperature. 



RESULTS AND DISCUSSION 

Figure ^a) shows the reduced longitudinal thermal 
conductivity k{T, B)/ k{T,0) as a function of the field 
applied parallel to the c— axis at different constant tem- 
peratures between 5 K to 20 K. In Fig.Q^b) we show the 
Hall resistance of the same sample. The decrease of k 
with field can be related to the decrease of the electronic 
contribution. The clear oscillations in k{B) observed at 
i3 > 1 T are apparently due to the quantization of the 
Landau levels and the crossing of the Fermi energy as 
the oscillations in the Hall effect indicate, see Fig. Q^b). 
Figure 13 shows the same data as Fig. ^ but in a linear 
field scale. 

It is noticeable the appearance of the plateau-like fea- 
tures at ~ 2 T and 4 T in the Hall resistivity that 
clearly suggests the occurrence of the Quantum Hall ef- 
fect (QHE) in graphite. In fact, the temperature depen- 
dence of the maximum slope (d|i?//|/di?)max vs. 
between two plateaus at 3.5 T, and measured to 70 mK 
shows a temperature dependence T"*^ with an exponent 
k = 0.42 similar to that found in QHE systems This 
result is not unexpected taking into account the quasi 2D 
structure of the sample. We note that the occurrence of 
the integral QHE in graphene has been predicted recently 
[2^ . The reason why it was not found before is related to 
the sample quality which affects the 2D behavior of the 
transport properties. Our studies show that the dimen- 
sionality is strongly affected by internal lattice defects, 
some of them appear to short circuit the graphene planes. 
The fact that we want to stress is the good correspon- 
dence between the features measured in k and Rh at the 
QL. To recognize this we show in the inset of Fig. ^b) 
the temperature dependence of the field at the onset of 
the plateau in Rh at ~ 3.7 T (o) and the position of the 
minimum in k (■). 

Figure 121 shows in more detail the field dependence of 
the basal-plane and Hall resistances, taken at 4.2 K, as 
well as of the thermal conductivity at 5 K. The ther- 
mal conductivity data shown in this figure were taken 
increasing and decreasing field. As seen in the figure, no 
significant hysteresis is observed. This result is in con- 
trast to the hysteresis observed by Ayache in his Ph.D. 
work ^17]. We speculate that a small temperature drift 
might have been the reason for the observed hysteresis. 

Can we understand the decrease with field and the os- 
cillations observed in k{B) within the Wiedemann-Franz 
relation? To answer this question we proceed as follows. 
We assume that the thermal transport of graphite is given 
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FIG. 1: (a) Reduced thermal conductivity Hi{B)/hi{0) vs. ap- 
plied magnetic field at different constant temperatures. The 
lines (2) and (3) were calculated from the Wiedemann-Franz 
law Eq. and Eq. Q at 10 K and 6 K using the mea- 
sured field dependent electrical resistivity for the same sam- 
ple. Curve (1) was obtained with the measured resistivity at 
10 K but with a lower Lorenz number L = 2.0 x 10"* WQK'^ . 
(b) Absolute value of the Hall resistance as a function of ap- 
plied field. The inset shows the temperature dependence of 
the field at the onset of the plateau in Rh at ~ 3.7 T (o) and 
the position of the minimum in k (■). 

by two contributions 0: 

K = Kp{T) + K,{T, B) , (1) 

where Kp is due to the phonons, the contribution of the 
atomic lattice with the appropriate lattice anisotropy, 
and Ke due to conduction electrons. Usually one assumes 
that the field dependence of the thermal conductivity is 
given only by the electronic part KeiT, B), which can be 
estimated with the WF relation. This universal relation 
relates the electrical resistivity p{T, B) with the thermal 
conductivity due to electrons by 
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FIG. 2: The same data as in Fig. 1 but plotted in a linear 
field scale. 
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FIG. 3: Longitudinal and Hall resistances (left scale) and 
the reduced thermal conductivity (right scale) as a function 
of magnetic field applied parallel to the c— axis. The two 
measured lines of the thermal conductivity data were taken 
increasing and decreasing field. Within the error there are no 
significant hysteresis. 



through the universal constant Lq = 2.45 x 10~® WfiK"^ 
at low enough temperatures. One can recognize easily 
the difficulty to measure accurately the field dependence 
of the electronic contribution to the thermal transport 
above ~ 4 K due to the small electronic contribution. 
From Eq. Q we expect for well ordered HOPG sam- 
ples at zero field a ratio between the electronic and total 
thermal conductivity Ke/n < 0.15 at 5 K, and < 0.05 at 
10 K. Literature data are in general in good agreement 
with these estimates ^3 01 • 

The relation (O holds strictly for elastic or quasielastic 
electron scattering and therefore the range of validity is 
usually set, either at low enough temperatures where the 
resistivity is temperature independent (impurity scatter- 
ing dominates), or at high enough temperatures where 
the electron-phonon scattering is large jl^- For the 
sample measured in this work, the temperature depen- 
dence of the electrical resistivity indicates a saturation 
below 10 K (curves for a similar sample can be seen in 
Refs. dill) and therefore at T < 10 K we expect to be 
roughly in the validity range of the WF-law. From the 
measured field dependence of the electrical resistivity we 
can calculate the relative change of the total thermal con- 
ductivity at a fixed temperature as: 



k{Q) 



k(0) 



1, 



(3) 



assuming that the phonon conductivity does not depend 
on magnetic field. In Fig. Q^a) we show three curves 
calculated with Eqs. (O and ^ using the measured p{B) 
at 10 K (curves (1) and (2)) and 6 K (3). Curve (1) 
was obtained with the same parameters as (2) but with 
Lq — 2.0 X 10~^WOK~^, assuming a decrease of Lq due 
to the possible influence of the inelastic scattering. 

From the comparison between the computed curves 
and the experimental ones one would tend to conclude 
that Eqs. (2) and (3) provide reasonably well the over- 
all decrease of k with field, within the geometrical errors 
in the measurement of both conductivities. Neverthe- 
less we should note that the WF law and Eq. © do 
not reproduce accurately the measured field dependence, 
see Fig. 1(a). Due to the electrical resistivity increase 
(a factor ~ 100 from zero field to 1 T at T < 10 K) 
the electronic contribution, according to the WF law, be- 
comes negligible. Therefore the electronic contribution to 
K should be negligible, e.g. at B > 0.4 T and T 10 K, 
and a saturation of k at larger fields is expected. This 
is not observed experimentally. This means also that 
the relatively small oscillations observed in the longitu- 
dinal resistivity above 1 T that accompany the features 
in the Hall effect (see Fig. 2) should not affect k accord- 
ing to Eq. ^ in contrast to the experimental results, see 
Fig. 1(a). These results clearly indicate that the WF law 
in its original form fails to explain the field dependence 
of the thermal conductivity in the QL of graphite and in 
the measured temperature range. 
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FIG. 4: The thermal conductivity difference between the min- 
imum (at B ~ 3.7 T) and maximum (at B ~ 5.5 T) Ak (see 
Fig.0 as a function of temperature. The right scale provides 
the estimate values per graphene layer (within a geometrical 
error of ~ 30%). 



For any future comparison with theory it may be use- 
ful to provide the absolute value of the oscillation am- 
plitude between minimum and maximum as defined in 
Fig. ^a). Figure 01 shows the experimental Ak obtained 
from Fig.^a) or Fig.|2Ia). It shows a tendency to satura- 
tion around Ak ~ 1.5 W/mK for T < 10 K. The oscilla- 
tion amplitude per graphene layer can be estimated from 
Ak^ ^ Ana/N, where a is the distance between layers 
and N the number of layers measured in parallel in our 
experiment within the 60 fim thickness of the sample. Us- 
ing this approximation we obtain Akl ~ 3 x 10~^^ W/K 
for T < 10 K. 

In the following we discuss our results taking into ac- 
count relevant work. A recently published theoretical 
work calculated the thermal conductivity of a 2D elec- 
tron gas at low temperatures and in a quantizing mag- 
netic field [i^ . Although some of the assumptions done 
there may not be valid for graphite, it is interesting to 
note that this work shows that the WF-law is violated for 
small Landau level broadening and at low enough tem- 
peratures when the difi^usion mechanism dominates. Ac- 
cording to this work the deviations from the WF law are 
due to the energy derivatives of the longitudinal electri- 
cal conductivity. At low enough temperatures and small 
level broadening the numerical results show a "two-peak" 
behavior, i.e. k as a function of magnetic field shows 
two maxima in the field region between two neighboring 
plateaus in the Hall resistance. A similar result has been 
obtained in previous theoretical work [sol IsH - In our 
case, however, we obtain the striking result that He{B) 
decreases in this field region whereas increases with field 
in the plateau region and reaches a maximum at the end 
of the corresponding plateau, see Fig. 01 If we would have 
localized QP in the field region of the plateau we would 
naively expect a decrease of Kg. On the other hand, the 



opposite behavior may be also possible, i.e. an increase 
of Ke with field, if the density of interacting QP would de- 
crease in this field region and the main scattering mech- 
anism is given by QP-QP scattering, as in the case of 
high-temperature superconductors [33|- In this case the 
theoretical description of the thermal conductivity may 
become more complicated to handle. 

According to recently pubhshed theory 0, H^, a 
magnetic field applied perpendicular to the graphene 
planes opens an insulating gap in the spectrum of Dirac 
fermions, associated with an electron-hole pairing, lead- 
ing to the excitonic insulator state below a field depen- 
dent transition temperature. The experimental value of 
the critical field of the field-driven metal-insulator tran- 
sition in graphite is ^^50 times smaller than the predicted 
in Refs. [23 ■ The discrepancy can be understood, how- 
ever, assuming that the Coulomb coupling, given by the 
dimensionless parameter g = l-ne^ jt^v (eg is the dielec- 
tric constant) 0, ' drives the system very close to the 
excitonic instability. In this case, the threshold field 
can be well below the estimated value of 2.5 T. The above 
analysis, together with the experimental evidence that 
only the perpendicular component of the applied field 
drives the MIT 3 1 appear to support the theoretical ex- 
pectations of a field-induced excitonic insulator state in 
graphite. In this case and according to Ref. [l^l we would 
expect a monotonic decrease of Ke(-B) with field and a 
kink at i3 ~ Be with a plateau region in the insulating- 
like state of the QP at B > Be- In the temperature range 
of our experiments the results do not show a clear kink 
at i? 0.1 T nor a plateau at higher fields and T > 6 K, 
see Fig. ^a), although one may tend to recognize it at 
T = 6 K and B ~ 0.5 T when the data are plotted in a 
linear field scale, see Fig.|5Ja). At the temperature limit 
of our system (5 K) the density of points at B < 1 T 
is too low to assure the existence or non-existence of a 
kink. Measurements at lower temperatures are needed 
to enhance the relative contribution of the QP to k and 
check whether the predicted feature is observable. This 
issue will be studied in the future. 

Can the oscillations in k{B > 1 T) be due to the lat- 
tice contribution Kp(B) via electron-phonon interaction 
as, for example, in antimony l35l |? The temperature de- 
pendence of K oc T^'^ at 5 K < T < 30 K speaks for 
phonon scattering by grain boundaries and not by elec- 
trons 0. The inelasticity parameter 77 = v/XlOc (here v 
is the sound velocity, A the magnetic length and tUc the 
cyclotron frequency), that provides an estimate of the ef- 
ficiency of the electron-phonon scattering, is ^ 0.01 at 
4 T for graphite. Therefore, unless there is an intersec- 
tion of Landau levels that favors acoustic phonon tran- 
sitions 36], it does not seem that the phonon-electron 
scattering can be significantly enhanced at high fields in 
graphite. The overall correlation of k(T, B) with the Hall 
resistivity indicates that the origin of the oscillations in 
k{B) should be related to a pure QP phenomenon. 
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In summary, high-resolution measurements of the mag- 
netic field dependence of the thermal conductivity in a 
quasi 2D sample of graphite show clear oscillations in 
the quantum limit. The Hall effect for the same sample 
shows quantum Hall effect features which are correlated 
to the features observed in k. With the measured lon- 
gitudinal electrical resistivity we show that the observed 
oscillations in n cannot be explained with the original 
WF law. Lower temperature measurements as well as an 
appropriate theoretical framework for graphite are highly 
desirable. 
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the thermocouple, M. A. Vozmediano and F. Guinea for 
fruitful discussions and their interest on this work, and E. 
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theoretical work. This research was partially supported 
by the DFG under DFG ES 86/6-3, FAFESF and CNFq. 
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